Comment on "Ferroelectrically Induced Weak Ferromagnetism by Design" 

The question of how ferroelectric polarization is coupled to magnetism in magnetoelectric 
multiferroics is of considerable current interest. A recent Letter [lj analyzes the important 
"ABCV class of perovskite multiferroics. A symmetry argument is presented that materials 
with antiferromagnetism on the B site lack linear magnetoelectric coupling Eplm = P • (L x 
M) in the free energy unlike A-site antiferromagnets [U [2]. Here P is polarization and L 
and M are antiferromagnetic and ferromagnetic moments. 

This Comment presents a distinct analysis of E PLM in ABO3 multiferroics. We show 
that the argument of Ref. |T| does forbid Eplm if the final low-symmetry phase contains only 
one distortion that, like P, breaks all inversion symmetries. In reality there are multiple 
distortions in this symmetry class, and cross-terms generate Eplm even in B-site materials, 
although the mechanism is different than in A-site materials. Additional differences can 
emerge between A-site and B-site materials under dynamical assumptions beyond the static 
free energy. 

For concreteness we use BiFe03, whose electronic structure has been extensively stud- 
ied j3jH]. Below HOOK it shows ferroelectric order with space group R3c and a 10-atom unit 
cell. In Ref. [TJ the free energy is expanded from a model R3c phase with inversion center 
on Fe. For now, assume that the polarization P arises from relative motion of the ionic 
sublattices [3j along the threefold axis z. There is a 13.8° rotation a of the O octahedra 
around Fe atoms (az is D in Ref. []]), which is even under Fe-site inversion I B and odd 
under Bi-site inversion I a- However, distortions beyond P and a must be present, either 
by counting R3c degrees of freedom or because the 1.4° counter-rotation of upper and lower 
triangles within an octahedron [5j [6] cannot be obtained by combining P and a. 

Coordinates for an oxygen distortion (3 (Fig. 1) that causes this counter-rotation and is 
orthogonal to P are provided in EPAPS. The (3 distortion, like P, is odd under I a and Ip- 
For A-site magnetism, we have inversion eigenvalues IaL = +L, and IpL = —L. For B-site 
magnetism, IaL = —L and IpL = +L. Allowed L x M terms are, with constants Ai, A 2 , 

F A = XiaP • (L x M) + X 2 a(3z • (L x M), 

F B = X x apP -(LxM) + X 2 az • (L x M). (1) 
More generally, P is some I A = I B = —1 distortion, and E PLM results from cross-terms 
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FIG. 1: Distortions a and (3 that reduce the ideal perovskite symmetry to R3c (along threefold 
axis). The distortion (3 combines with a at order af3 to give both the observed counter-rotation 
and a linear magnetoelectric coupling. 

between P and other such distortions, here (3. Magnetoelectric coupling should be sought 
in B-site structures with large /?-type distortions as well as in A-site structures with large 
a. 
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EXPLICIT CONSTRUCTION OF a AND (3 DISTORTIONS (APPENDIX TO COM- 
MENT ON "FERROELECTRIC ALLY INDUCED WEAK FERROMAGNETISM 
BY DESIGN") 

It is convenient to use the following notation: Vectors written without parenthesis denote 
components along the primitive vectors: x,y,z = xa + yb+ zc. A convenient set of primitive 
vectors for the 10-atom R3c unit cell is a = (0, 1, 1), b = (1, 0, 1), and c = (1, 1, 0). Here 
vectors with parenthesis are in a cartesian basis, and for convenience we choose the origin 
to be on top of a B atom. Note that if Z, to, n is a set of integers, our notation implies that 
x + l,y + to, z + n is equivalent to x, y, z apart from a Bravais vector translation. 

For a general AB0 3 perovskite with R3c symmetry, both A and B atoms are located at 
Wyckoff 2a positions: The A atoms are located at x 1 , x', x' and x' + 1/2, x' + 1/2, x f + 1/2, 
while the B atoms are located at x" , x" ', x" and x" + 1/2, x" + 1/2, x" + 1/2. 

The six oxygens must be located at Wyckoff 6b positions: 



We may parametrize a set of paraelectric R3c structures (under the definition of polar- 
ization assumed for this example in the text, in which the polarization is well described as 
translation of the anion and cation lattices; if the definition of P changes, then that of (3 
changes, but by counting such a (3 always exists) with two degrees of freedom a, (3 according 
to: 



x,y,z 



z,x,y y,z,x 
ill 111 





This leads to the following atomic positions: 



Ra 



1 1 1 

4'4'4 



3 3 3 
4' 4' 4 
1 1 1 

2' 2' 2 



Rb 




Rq 



— - + a + f3, - — a + (3, - — 2/3 and cyclic permutations 
5 13 

a + (3, — h<^ + /3, 2/3 and cyclic permutations. 
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Note that any equivalent lattice site could equally well be chosen for the oxygen locations, 
e.g., the choice above — |, |, \ is equivalent to— |, — |, |, or |, |, |, and so on. 

The a and /? distortions defined in this way are non-polar, i.e., orthogonal to P. To see 
this, consider the electric dipole moment produced by the oxygens, 

Po = Qo^2Ro= ^[3 + 2(/ + m + n)] x 1,1,1. (2) 

O l,m,n 

This is seen to be independent of a, (3 and is equal to zero when summed over a symmetric 
set of lattice vectors (the electric dipole moment of the A and B atoms is also zero after a 
suitable summation). Direct computation shows that there is a counter-rotation by an angle 
proportional to a/3, as shown in Fig. 1 of the Comment. 

The transformation properties under inversion on top of an A atom, 7^, and inversion on 
top of a B atom, I B) are straightforward to compute from the atomic positions above. For 
Ib ) just apply inversion directly to the .Ro coordinates, while for I a we must first subtract 
1/4, 1/4, 1/4, invert, and then add 1/4, 1/4, 1/4. The final result is simple to state: When 
I a is applied, the original structure is recovered if we apply the transformations a — > — a 
and (3 — > —f3. Under I B , the original lattice will be recovered only if we apply a —> +a, and 
(3 —> —f3. In summary, 

I A (aJ) = (-a,-p) 

I B (aJ) = (a, -(3). (3) 
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